Resonances in rotationally inelastic scattering of OH(X 2 U) with helium and neon 
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We present detailed calculations on resonances in rotationally and spin-orbit inelastic scattering 
of OH (X 2 II, j = 3/2, Fi, f) radicals with He and Ne atoms. We calculate new ab initio potential 
energy surfaces for OH-He, and the cross sections derived from these surfaces compare favorably 
with the recent crossed beam scattering experiment of Kirste et al. [Phys. Rev. A 82, 042717 
(2010)]. We identify both shape and Feshbach resonances in the integral and differential state-to- 
state scattering cross sections, and we discuss the prospects for experimentally observing scattering 
resonances using Stark decelerated beams of OH radicals. 



I. INTRODUCTION 

Measurements of state-to-state cross sections pro- 
vide important tests of the reliability of computed po- 
tential energy surfaces (PES's) describing the interac- 
tion of atoms and molecules Cross sections for 
collision-induced rotational transitions are sensitive to 
the anisotropy of the PES. Since non-bonding interac- 
tions are relatively weak, the magnitudes of the cross 
sections are mostly sensitive to the repulsive part of the 
PES, except at very low collision energies. An alterna- 
tive, spectroscopic approach to gaining information on 
PES's is the determination of the energies of the bound 
levels of van der Waals complexes of the collision part- 
ners @, 0] . The energies of the bound levels are mainly 
sensitive to the attractive part of the PES's. As we go 
up higher in the manifold of these weakly bound lev- 
els, the energies of these levels eventually become higher 
than the dissociation energy of the complex, and such 
levels are quasi-bound. These quasi-bound levels are of- 
ten described as resonances and can be thought of as a 
distortion of the continuum in the collision energy de- 
pendence of state-to-state cross sections [||. In inelastic 
scattering, resonances are called shape or orbiting res- 
onances when the quasi-bound levels involve monomer 
levels that are the same as in the initial or final level of 
the collision-induced transition, or Feshbach resonances 
when they involve different monomer states 0, H| . Due 
to their sensitivity to the PES, resonances can reveal im- 
portant information on the PES [HQ. So far, however, 
resonant structures in scattering cross sections have been 
experimentally observed only in exceptional cases 

The crossed molecular beam technique has been an 
extremely useful tool for the determination of state-to- 



state cross sections, both integral and differential, as well 



as their dependence upon the collision energy The 
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recently developed Stark deceleration technique, taking 
advantage of the interaction of polar molecules with time- 
varying electric fields, has allowed continuous tuning of 
the beam velocity [l2| ■ This has facilitated measurements 
of the collision energy dependence of state-to-state in- 
tegral cross sections down to energies of 70 cm -1 [l3| . 
Moreover, the velocity spread in such decelerated beams 
is much smaller than in conventional molecular beams. 
Thus far, an energy resolution of > 13 cm -1 has been 
achieved for collisions of OH radicals with rare gas atoms 
[TH-EH . This resolution is mainly limited by the velocity 
and angular spread of the atomic collision partner, and is 
too low to experimentally resolve scattering resonances. 
A recent study has shown that the energy resolution can 
be improved significantly by an appropriate choice of the 
beam velocities and interaction angle [l^ . When these 
measures are put into practice in the laboratory, collision 
energy resolutions can be obtained that may enable the 
observation of scattering resonances. 

Atom-molecule collisions are the simplest type of col- 
lision process in which rotationally inelastic transitions 
can be observed. Early calculations [l?], [HI on rotation- 
ally inelastic scattering of N2 molecules with He atoms 
have shown that resonances occur at low collision ener- 
gies, but the experimental verification of these predic- 
tions was not yet possible. Collisions of OH(A 2 n) with 
rare gases have emerged as paradigms of scatte ring of an 
open-shell molecule with an atom [H, [l5, 19 24]. The 
OH-rare gas systems are good candidates for the obser- 
vation and analysis of resonances in rotationally inelastic 
collisions because the collision energy can be reduced by 
Stark deceleration of the OH beam. Since OH(A 2 n) is 
an open-shell molecule with orbital degeneracy, the col- 
lision dynamics is governed by two PES's, [25[ and in- 
teresting multi-state dynamics can occur. Of particular 
interest for the study of resonances are the OH-He and 
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OH-Nc systems, since the dissociation energies of these 
systems are smaller than the rotational level spacings of 
the OH radical. The resonance features associated with 
the various rotational levels are therefore well separated. 
The shallow van der Waals wells support only one or two 
stretch vibrational levels [26j], resulting in a rather sim- 
ple, yet interesting, analysis of the resonances. Shape 
res onances in OH-He collision s were previously analyzed 
by IDagdigian and Alexandeilp i | in a study of elastic de- 
polarization. Bound states of the OH-He complex have 
been investigated spectroscopically by Han and Heaven, 
who identified complex features as scattering resonances 
in OH(^)-He H3- 

Here, we present a detailed and precise study of scat- 
tering resonances in the OH-He and OH-Ne systems, in 
order to develop insight into the nature and strength of 
the resonances and to assist in the experimental search 
for such scattering resonances. For the OH-He system, 
we have calculated new three dimensional potential en- 
ergy surfaces. The collision energy dependence of the 
relative state-to-state integral scattering cross sections 
that are derived from these potentials compares more 
favorably with recent experiments [l5j than the results 
from previous calculations. For the correct assessment 
of the resonances, the calculations are performed on a 
very fine grid of collision energies, and particular care is 
taken to converge the calculations to avoid numerical ar- 
tifacts to be interpreted as resonant structures [28j]. We 
characterize the resonances with various techniques, in- 
cluding the adiabatic bender model [H, [30| and collision 
lifetime analysis [3l|. We investigate how the differen- 
tial cross section for several transitions changes as the 
collision energy is scanned through the resonances, and 
observe dramatic effects. 

This paper is organized as follows: The details of the 
scattering calculations are briefly presented in Section [TT1 
In Section lHll wc describe the new three-dimensional (3D) 
PES's that are developed for OH-He. Section [TV] de- 
scribes our calculations on the state-to-state scattering 
cross sections in OH-He collisions. A detailed analysis 
of shape and Feshbach resonances is given. Section [V] 
presents similar results for the OH-Ne system. A dis- 
cussion of the prospects for observing these resonances 
in crossed beam experiments using a Stark decelerator, 
by either recording the integral or the differential cross 
sections, follows in Sec. IVII 



II. SCATTERING CALCULATIONS 

The theory of scattering between a molecule in a 2 n 
electronic state and a structureless atom is well estab- 
lished [25| . The interaction can be described by two 
PES's corresponding to states of A' and A" symmetries. 
For OH-He, we have constructed new PES's, which are 
explained in Sect. IIIIl wh ile for OH-Ne, we used the 
PES's by ISumivosh i et al.l [32l ] . Close-coupling calcula- 
tions were performed both with the HIBRIDON suite of 
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FIG. 1: Energies of the lower rotational levels of OH(X 2 II). 
The A-doublet splitting is exaggerated for clarity. The initial 
level for all scattering calculations is the j = 3/2, Fi, f level. 



programs [33| , and with a second independent scattering 
program for open-shell diatom-atom scattering described 
in Ref. [HI- Care was taken to independently check the 
results with the two scattering programs and to converge 
the cross sections. For OH-He the maximum total an- 
gular momentum was J = 100.5 - 140.5, depending on 
the collision energy, and the channel basis consisted of 
all rotational levels of OH with j < 6.5, while for OH- 
Ne the channel basis consisted of all rotational levels 
with j < 7.5. In this paper, we calculate cross sections 
from fully converged close-coupling calculations in order 
to study resonances in inelastic collisions between low- 
lying rotational states of the OH radical. For reference, 
the energies of the lower rotational levels of OH(A 2 n) 
arc displayed graphically in Fig. [T] 



III. 3D OH-HELIUM POTENTIAL 

A crucial role in the scattering calculations is played 
by the interaction potential. In Ref. HH a detailed ex- 
perimental and theoretical study of inelastic scattering 
between OH radicals and the rare gas atoms He, Ne, Ar, 
Kr and Xe was performed. The theoretical results in 
that study were shown to be in excellent agreement with 
experimentally measured inelastic cross sections. The 
agreement between theory and experiment was, although 
still very good, the worst for the OH-He system. It was 
believed that this was due to the quality of the PES, since 
for the OH-He system a smaller basis set was used in the 
calculations than for the other systems. For this reason, 
we construct here a new potential for the OH-He system. 
We note that in Ref. [34| the experimental resolution was 
unfortunately not yet high enough to observe resonances. 

In trying to improve the agreement with the experi- 
mental results, we first constructed new 2D PES's for 
OH-He collisions. This was done by enhancing the ba- 
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FIG. 2: Illustration of the OH molecule and the He atom 
containing the relevant coordinates as used in the fitting of 
the OH-He potential. R is the length of the vector R that 
connects the He atom and the center-of-mass (Q) of the OH 
molecule, while 9 is the angle between R and the OH bond 
vector r (pointing from O to H) of length r. R a is the length 
of the vector R a that connects the He atom and the H atom, 
while 9 a is the angle between R a and r. Rt is the length of 
the vector R;, that connects the He atom and the O atom, 
while 9i, is the angle between R5 and — r. The point X marks 
the location of the midbond orbitals. 

sis set for the coupled-cluster calculations of the interac- 
tion energy from the augmented tripl e-zeta corr elation- 
consistent basis set (AVTZ) used by iLee et all [26[ to 
the quintuple-zeta basis set (AV5Z). We computed the 
interaction energies with the open-shell single and dou- 
ble excitation coupled cluster method with perturbative 
triples as implemented in the MOLPRO package [HI . The 
interaction energies were evaluated for 288 geometries on 
a two-dimensional grid with 12 Gauss-Legendre points 
in the Jacobi-anglc 9. The OH bond length was fixed 
at the v ibrationall y averaged distance of ?'o = 1-8502 ao, 
whereas ILee et ail used the equilibrium distance r e . The 
relevant geometry is illustrated in Fig. O We included 
midbond orbitals (3s.3p,2d,lf.lg) with the exponents of 
Ref. These midbond functions were centered on the 
vector R that connects the He atom and the center-of- 
mass of the OH molecule, at a distance from the helium 
atom that is half the distance of the helium atom to the 
nearest atom of the OH molecule. Also the counterpoise 
correction of Boys and Bernardi was applied [13 ■ The 
grid of atom-molecule separations consisted of 18 points 
ranging from R = 3 ao to 9 ao at short range and 6 points 
on an approximately logarithmic scale up to 25 ao at long 
range. 

As mentioned before, two potential energy surfaces be- 
longing to states of A' and A" symmetry are involved in 
the OH-rare-gas atom scattering. The average and 
half-difference Vd of these potentials can be expanded in 
Racah normalized spherical harmonics Ci >mi namely 

^ = Va> + V*' =f\ i|0 (fl)C t , o (g,0), 
1=0 

Vd = Va "~ Va ' =^v h2 (R)C l>2 (9,0), (1) 
1=1 

where we included all terms up to l max = 11. In the long 
range (R > 10ao), the expansion coefficients vi o(R) were 



fitted to inverse powers R " with n > 6, namely 

v* (R,e)= £) Q, n Mp, (2) 

71=710 (0 

where we note that the allowed values for n depend on 
I [H}. For example, for I = we have no(l) = 6 and 
only even values of n are present, while for I = 1 we have 
no(l) = 7 and only odd values of n are present. From 
the fitted coefficients we only kept the leading long-range 
terms for I = to / = 4. We used the Tang-Tocnnics 
damping function [39l ] 

/»(*) = 1-«^Ef ^ 

to damp these five long-range terms in the short range 
with j3 = O.GaQ 1 . In the short range (R < 5.5ao), the ex- 
pansion coefficients v^o(R) were fitted to an exponential, 
namely 

vf(R) = Sl e~ a ' R . (4) 

The difference between the ab initio interaction energies 
and the analytic long range and short range functions was 
fitted with a reproducing kernel Hilbert space (RKHS) 
method The RKHS parameter m was chosen such 
that the RKHS fit would decay faster than the leading 
long-range term for each I. The RKHS smoothness pa- 
rameter was set to 2. For the expansion coefficients of 
the difference potential, vi^(R), no analytic short range 
and long range fit was performed, so that everywhere 
the RKHS method was used. Using the described proce- 
dure, we obtained an accurate fit to the ab initio points. 
More details of the fit can be found on EPAPS [4l[ , where 
we provide a FORTRAN 77 code for the two-dimensional 
AV5Z potential. 

We found that the absolute minimum of the fitted po- 
tential is located at 6 = 68.7°, R = 5.69 a on the A' PES, 
corresponding to an interaction energy of Va> = —29.8 
cm -1 . The minimum potential energy values for 9 = 0° 
and 9 = 180° were found at R — 6.56 ao and R = 6.09 ao, 
giving rise to Va'/a" = —27.1 cm -1 and Va> ia" = —21.6 
cm , respectively. Fo r compari s on, we also mention 
the values obtained by ILee et all [261 ] , who found that 
the absolute minimum of their potential was located at 
9 = 68.6°, R = 5.69 ao for A' symmetry, with an interac- 
tion energy of Va> = —30.0 cm . The minimum values 
for 9 = 0° and 9 = 180° were found at R = 6.54 a 
and R = 6.09, giving rise to Va>/a" = —27.1 cm -1 and 
Va 1 /A" = —21.8 cm -1 , respectively. The two potentials 
are seen to give very similar results for the local and 
global minima. Moreover, using the new AV5Z potential 
for scattering calculations, we found only a very slight im- 
provement in the agreement with the experimental data. 

Therefore, we tried to improve the PES further by tak- 
ing the vibrational motion of the OH radical into ac- 
count. To this end, we computed the interaction ener- 
gies of the OH-He system on a three-dimensional grid. 
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At short and intermediate range we used a step size of 
AR = 0.25 a for 3a < R < 12.5 a and Ar = 0.25 a 
for 0.75 ao < r < 4.5 ao- For the angle 9 we used an 
equidistant grid of 16 points including and 180° with a 
spacing of AO = 12°. At long range we used 4 equidistant 
points between 14 < R < 20, while we used a step size 
of Ar = 0.5 ao for 0.75 ao < r < 4.25 ao. The distance 
r = 4.5 ao was also included in the long-range fit. For 
the angle 9 we used an equidistant grid of 9 points with 
a spacing of A9 = 22.5°. On this grid we computed the 
interaction energies with a triple-zeta basis set (AVTZ) 
and usin g m idbond orbitals with geometry-dependent ex- 
ponents [42|. Especially for large r and small R, the 
electronic structure calculations did not always converge. 
Then, we obtained the energy for the corresponding grid 
point by means of interpolation or extrapolation from 
neighbouring grid points. 

To perform the fit of the sum interaction potential V s , 
we proceed in the following way. We represent the po- 
tential as a sum of three terms, namely 

V s (R,e,r) = 

V s sr (R a ,9 a ,r) + V s SI (R b , 9 b , r) + V s h (R, 9, r), (5) 

where the different coordinates are defined in Fig. [2] 
This representation is convenient because the coordinates 
of the first and second term of Eq. ([5|) are ideally suited 
to describe the short-range behavior near the H and O 
atom, respectively, while the coordinates of the third 
term are convenient to describe the long-range behavior. 
The short-range terms arc fitted by 

'max 

^ sr (i?^,r) = Y j e-^ Ri P l {cos9 i )sf ) 

1=0 

'max ^max n max 

+ E E E Rfe-^Piicose^e-^s^ (6) 

(=0 k=0 n=0 

where i = a, 6, while Pi(x) are Legendre polynomials 
corresponding to the functions Cifl(0,0) of Eq. [1] We 
used the values i max = 1, n^ ax = 3, fc° iax = 8, ^ iax = 7, 
n max = 3, fcfnax = 8 and ^f nax = 5. The long range term 
is fitted by 

V>W,r) = f^ 7 J2il^l Pl{c0S 9)c nl (r\ (7) 

n=6 1=0 

where /„ is the damping function of Eq. ([3]). Nonzero 
values of c„; occur only for even values of I + n, and then 
they are given by 

3 

c nl (r)=c° nl +J2r k e- a - r3 c nlk . (8) 

We use two different values for a n , namely a„ = ol\ for 
6 < n < 9, and a n — an for 10 < n < 13. For the dif- 
ference potential Vd similar fit functions are used, only 
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FIG. 3: A' potential energy surface. The OH radical lies on 
the horizontal axis, with the center-of-mass of the molecule at 
the origin. The O atom lies left of the origin, the H atom to 
the right. For each geometry of the complex, defined by the 
OH bond length r and the position (ih b , J/He) of the He atom, 
the interaction energy is calculated, resulting in contours with 
the unit of cm -1 . The three plots differ in the OH bond 
length, namely in panel a) we have r — 0.75 ao, in panel b) 
r = 1.8324 ao, and in panel c) r = 3.00 ao. 



now the Legendre polynomials Pi(x) are replaced by as- 
sociated Legendre functions Pf{x) corresponding to the 
Racah spherical harmonics (7(^(6*, 0) of Eq. [TJ so that 
all sums start with I = 2. Moreover, we use /^ ax = 2, 
,l max = 5, fc max 5, t max = 6, rt max = 5, fc max = 4 and 
ljjj ax = 6. The linear and the nonlinear fit parameters 
were determined by minimizing a weighted least-squares 
error. 

By evaluating the analytic representation of the poten- 
tial on the ah initio grid, we were able to compare the 
fitted energy values with the ab initio values. We found 
that we only obtained a reliable fit for OH bond lengths 
r < 3dQ. At smaller values of R the largest relative error 
of an analytic value compared to an ab initio value for 
r < 3oq was 6.67% for the sum potential and 1.10% for 
the difference potential. At large i?, again considering 
only r < 3ao, the largest relative error was 3.23 % for 
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FIG. 4: Similar to Fig. [3l but showing the A" potential energy 
surface. 



the sum potential and 3.32% for the difference potential. 
We also calculated the average relative error, which for 
the sum potential was 0.32% at short range and 0.49% 
at long range, while for the difference potential it was 
0.04% in the short range and 0.90% in the long range. 
In Fig. [3J we show two-dimensional contour plots of the 
fitted OH-He A' PES for r = 0.75 a , r = 1.8324 a and 
r = 3.00 Oo, while in Fig. [4] the same plots are shown 
for the fitted A" PES. For the equilibrium bond length 
r e = 1.8324 ao, our fit of the 3-dimcn sional pote ntial is in 
very close agreement with the PES of lLee et al.l (26j , as it 
should, since both PES's were calculated with the same 
ab initio method using the same basis set. The absolute 
potential energy minimum for r = 1.8324 ao is located at 
9 = 69.2°, R = 5.69 ao f° r A' symmetry, leading to an 
interaction energy of Va> = —30.0 cm -1 . The minimum 
values for 9 = 0° and 9 = 180° were found at R = 6.55 ao 
and R = 6.09, giving rise to Va'/a" — —27.2 cm' 1 and 
Va'/a" = -21.7 cm -1 , respectively. A Fortran 77 code 
to generate the interaction potential V(R,9,r) is made 
available as an EPAPS document [4l[ . 

To use the three dimensional potential for scattering, 
we started with the three-dimensional AVTZ potential, 
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FIG. 5: Relative state-to-state inelastic scattering cross sec- 
tions of OH (X 2 H,j = 3/2, Fi, f) radicals with He atoms. 
The experimental data points from Ref. [34| are shown as dots, 
while the theoretically calculated cross sections with the po- 
tential of lLee et alJ [2(| are included as solid curves, and the 
results with the adiabatic potential as dashed curves. On the 
vertical axes of the plots, 100 % refers to the total inelastic 
cross section. Relative cross sections for inelastic collisions 
populating the (a) j = 3/2, Fi, e (black), j = 5/2, Fi, e (pur- 
ple), and j = 5/2, Fi, / (orange) states; (b) the j = 1/2, Fb, e 
(brown) and j = 3/2, F2, f (pink) states; (c) the j = 1/2, Fi, f 
(red) and j = 3/2,F 2 ,e (blue) states; (d) the j = 7/2, F\,e 
(green) and j = 7/2, Fi, / (cyan) states. 



and then subtracted the values of this potential at r = ?*o 
and added the two-dimensional potential calculated for 
r = tq at the AV5Z level. This implies that the depen- 
dence of the intermolecular potential on the most rele- 
vant coordinates for the scattering calculation, namely 
R and 9, is computed at the AV5Z level for r = ro, while 
the variation of the potential with the OH bond length 
r is taken into account at the AVTZ level. We solved 
for the OH vibrational motion in the full 3D potential 
generated by taking the intermolecular potential energy 
V{R,9,r) and adding the free OH monomer potential 
Von(r) [43[. For fixed R and 9 this leads to an effec- 
tively one-dimensional problem that can be easily solved 
by standard numerical methods, such as the discrete vari- 
able representation based on sinc-functions (sinc-DVR) 
[44]] . Taking the resulting ground state energy for each 
R and 9 and subtracting the v = monomer vibrational 
energy in the absence of the He atom then results in an 
adiabatic two-dimensional PES. We found, actually, that 
this adiabatic potential is very similar to the 'diabatic' 
potential obtained by first calculating the lowest vibra- 
tional state of OH in the monomer potential Voh (?*) and 
then averaging the interaction potential V(R, 9, r) over 
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this ground state. The two methods are expected to give 
similar results since the vibrational levels of OH are well 
separated in energy, so that the weak OH-He interaction 
gives only a slight admixture of the higher vibrational 
states of OH. 

The inelastic OH-He cross sections with OH initially in 
the j = 3/2, Fx, f level were calculated with the adiabatic 
potential. The results are shown in Fig. [5j where the ex- 
perimental data of Ref. is also shown, as well as the 
scattering results obtained with the potential of Lee et al. 
pol ] . The theoretical data are convoluted with the exper- 
imental energy resolution. To this end a Gaussian energy 
distribution is taken with a standard deviation that is a 
function of the energy. The value of the standard devi- 
ation ranges from 24 cm -1 at low collision energies to 
59 cm -1 at the highest collision energies. We note that 
the relative cross sections are plotted, rather than the 
absolute cross sections, because these relative cross sec- 
tions are experimentally measured. More details can be 
found in Refs. HH, We see from Fig. [5] that the over- 
all the agreement with experimental data has improved 
noticeably with the adiabatic potential. 



IV. OH-HELIUM COLLISIONS 
A. State-to-state integral cross sections 

For OH-He collisions, the state-to-state scattering 
cross sections were calculated with the adiabatic poten- 
tial described in the previous section. In Fig. [SJ the en- 
ergy dependence of the state-to-state integral cross sec- 
tions for several transitions out of the j = 3/2, Fx, / level 
of OH are shown. This level, which is the higher A- 
doublet component of the ground rotational level (see 
Fig. [TJ, can be selected with the Stark decelerator since 
it is low-field seeking in an inhomogeneous electric field 
fl3| . The cross sections are computed on a very fine grid 
of energies to be able to study resonant features in detail. 

Away from the resonances, these results are in good 
agreement with those previously reported by Klos et al. 
[2ll ]. As noted by these authors, there is a propensity for 
transitions preserving the total parity. The cross sections 
are found to be smaller for transitions with large energy 
gaps. The initial and final levels of the two transitions 
shown in Fig. EJa) have a rather large energy separation 
(> 100 cm -1 ), and the total parity is inverted during the 
transitions. Hence, the cross sections for these transitions 
are small. 

Resonances can be observed in Fig. [5] near the colli- 
sion energies corresponding to thresholds for excitation 
of the OH radical to higher rotational and spin-orbit lev- 
els. Both shape resonances, which appear right above 
the threshold energies for the final levels, and Feshbach 
resonances, which appear near the energies where higher 
rotational levels than the considered outgoing channel be- 
come open, are present. Except for the j = 3/2, Fx, f — > 
j = 3/2,Fx,e and the j = 3/2, Pi,/ j = 1/2,F 2 ,/ 
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FIG. 6: State-to-state integral cross sections vs. collision en- 
ergy for transitions out of the OH j — 3/2, Fi, f level in colli- 
sions with He. The final levels are indicated for each transition 
for which the cross section is plotted. 



transitions, the Feshbach resonances arc not significant 
compared with the background continuum. The j = 
3/2, Fi, / — > j — 3/2, Fi, e transition dominates at low 
collision energies and also gives rise to shape resonances 
with cross sections peaking above 10 A 2 . However, these 
shape resonances occur at collision energies of only a few 
wavenumbers. 

In the following subsections, we analyze the shape res- 
onances in the j = 3/2, Fi, / — > j = 5/2, Fi, e transition 
and the Feshbach resonances in the j = 3/2, Ft, / — > j = 
1/2, F2, f transition. The former transition has a large 
cross section; the latter transition exhibits strong reso- 
nances that show the largest enhancement compared to 
the background. 
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FIG. 7: (a) State-to-state integral cross sections vs. collision 
energy for the j = 3/2, Fi, f —¥ j = 5/2, Fi, e transition of OH 
in collisions with He (thick line), and the integral cross section 
convoluted with Gaussian energy distributions of FWHM of 
1 and 5 cm -1 (thin lines), (b)-(e) Differential cross sections 
dcr/dD. of the above transition at several energies marked as 
dashed lines with Roman numerals in (a), together with the 
differential cross sections convoluted with Gaussian energy 
distributions as in (a). 



B. Shape resonances 

Shape resonances result from quasi-bound states of 
the van der Waals complex formed by the collision part- 
ners at energies just above the threshold for the final 
level. All integral cross sections plotted in Fig. [5] dis- 
play shape resonances. In this subsection we analyze the 
shape resonances associated with the j = 3/2, F\,f — > 
j = 5/2, F\, e transition since it has a large integral cross 
section. Figure [3a) displays these resonances on an ex- 
panded energy scale. Several maxima, with increasing 
peak width vs. energy, can be observed, as was also pre- 
viously found_p4| for OH-He and other He-molecule sys- 
tems |5l.l45l.l46l|. 

To gain more insight, we employ the adiabatic ben- 
der model [1^, [3(| to analyze the shape resonances. The 
method is similar to the previous analys is of OH-He col- 
lisions by iDagdigian and Alexander] [24j . except that we 
used a close-coupling channel basis instead of a coupled- 
states one. The full Hamiltonian with the inclusion of 
Coriolis coupling and only the radial nuclear kinetic en- 




FIG. 8: Plots of the OH-He adiabatic bender curves that cor- 
relate with the OH j = 5/2, Fi,e level, obtained from close- 
coupling calculations. Curves are labeled with where J, 
p, and n are the total angular momentum, the total parity of 
the scattering wavefunction, and the cardinal index, respec- 
tively. 



ergy excluded is diagonalized as a function of R. The 
eigenvalues define a set of adiabatic bender potential en- 
ergy curves, which are labeled by the total angular mo- 
mentum </, the total parity p of the scattering wavefunc- 
tion, and the cardinal index n of the eigenvalue. In this 
paper, we will use the symbol J^ + ' and to label 

close-coupling adiabatic bender curves with p = +1 and 
p = — 1, respectively. 

Figure [5] shows several adiabatic bender curves that 
correlate with the OH j = 5/2, Fx, e level. The curves 
marked with 5/2^ and 3/2^ are the two lowest ly- 
ing adiabatic bender curves, each of which supports only 
one bound stretch level, with energies of 77.47 cm -1 and 
78.15 cm -1 , respectively As J and n increase, the curves 
move up in energy and the well depths become smaller. 
As a consequence, some of the bound levels become quasi- 
bound, and for the high lying curves (for example, the 

curve shown in Fig. |SJ) the wells are too shallow 
to support any quasi-bound levels. 

To compute the energies of the shape resonances, we 
treat the adiabatic bender curves in conventional one- 
dimensional scattering problems and calculate the phase 
shift. We should be able to observe rapid changes by 
7r, signifying resonances, in the collision energy depen- 
dence of phase shift Q. Fig. |H] shows the phase shift as 
a function of collision energy for all the adiabatic bender 
curves that have such a feature. We see from Fig. [5] that 
resonances in six adiabatic bender curves contribute to 
each of the three peaks shown in Fig. [7{a) (labeled as I, 
II, and III), which occur at 84.8, 87.6 and 91 cm -1 , re- 
spectively. The resonance features in Fig.[7][a) and Fig. [3] 
match well both in energy and width. Note that in or- 
der to distinguish the phase shift in different adiabatic 
bender curves, Fig. [HJb) and (c) do not show the whole 
range of the resonances, and thus the resonant changes 
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FIG. 9: Phase shifts as a function of collision energy for OH(j = 5/2, Fi, e)-He collisions, obtained from close-coupling adiabatic 

bender curves described in the text and shown in Fig. [8] Curves are labeled with \ where J, p, and n are the total angular 
momentum, the total parity of the scattering wavefunction, and the cardinal index, respectively. 



in phase shift shown are less than tt. 

It is also interesting to compare the differential cross 
sections for collision energies on and off a resonance. Fig- 
ures [3Jb)— (e) display differential cross section for several 
energies marked in Fig. [TJa) with Roman numerals. The 
center of the two major peaks are marked as I and II, 
while III and IV correspond to non-resonant energies. 
We observe significant backward scattering for energies I 
and II, likely because of an increased time delay of colli- 
sion due to the formation and decay of quasi-bound levels 
of the van der Waals complex. Backward peaks are in- 
significant for energies III and IV. We will further discuss 
this topic in the next subsection. 



C. Feshbach resonances 

In Feshbach resonances, quasi-bound levels of the colli- 
sion complex associated with a given rotational level dis- 
sociate to yield the molecule in a lower-energy rotational 
level. We consider here Feshbach resonances associated 
with the j = 3/2, F x ,f ->• j = 1/2, F 2 , / transition. This 
transition was chosen for detailed study since the reso- 
nance features show a ca. 4-fold increase over the contin- 
uum background (see Fig.|6j). Figure ITOfa) displays these 
resonances on an expanded energy scale. 

It is seen that a rich set of Feshbach resonances exists 
in a collision energy range of several cm -1 below the ener- 
getic threshold for opening of the F-2, j = 3/2 level at 188 
cm -1 . Figure [TT1 displays the contribution to the integral 
cross section for the j = 3/2, Fi, f — > j = 1/2, F 2 , f tran- 
sition from various values of the total angular momentum 
J (partial cross sections). The individual partial cross 
sections exhibit one or more peaks, and their energies 
shift toward higher collision energy as J increases. For 
J > 13/2, no significant resonances can be found in the 
energy dependence of partial cross sections. 

We performed an adiabatic bender analysis similar to 
that described in subsection IIVBI We calculated adia- 
batic bender potentials by diagonalizing the Hamiltonian 
expressed in a close-coupling channel basis. Since all pos- 
sible values of I (the orbital angular momentum of the 
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FIG. 10: (a) State-to-state integral cross section vs. collision 
energy for the j — 3/2,Fi,f — > j — 1/2, Fa, / transition of 
OH in collisions with He (black line), and the integral cross 
section convoluted with Gaussian energy distributions with 
FWHM of 1 (red line) and 5 cm" 1 (blue line), (b)-(j) Differ- 
ential cross sections da/di} of the above transition at several 
energies marked as dashed lines and Roman numerals in (a), 
together with the convoluted differential cross sections with 
Gaussian energy distribution as in (a). 
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j = 3/2, Fi, f — s> j = 1/2, F2, f transition of OH in collisions 
with He for total angular momentum J < 11/2. The dotted 
vertical lines denote the computed energies of the van der 
Waals stretch levels supported by the close-coupling adiabatic 
bender curves. 



van dcr Waals complex) are included in the channel ba- 
sis, there are multiple adiabatic bender curves for each 
value of J. These adiabatic bender curves look similar 
to those shown in Fig. [5] and are not plotted here. The 
energies of the van der Waals stretch levels supported by 
these curves were derived using a fixed step-size discrete 
variable representation (DVR) method |47l |48|. To treat 
levels that might be quasi-bound, an infinite barrier was 
placed at the maximum of the centrifugal barrier on each 
adiabatic bender potential. For curves associated with 
large J, this approximation will lead to calculated ener- 
gies higher than they should be and could lead to signifi- 
cant error. These computed energies are shown as dotted 
lines in Fig. [IT] There is a reasonable match between the 
energies of the resonances and of the bend-stretch levels, 
especially for small J. 

Fig. [TUTbWi) display the differential cross section for 
the OH j = 3/2, F u f ->• j = 1/2, F 2 , f transition at sev- 
eral energies marked on Fig. llOf a) with Roman numerals. 
The energies at I and IX are not at a resonance, and the 
differential cross sections show little backward scattering, 
while II - IV correspond to resonance energies, for which 
some backward scattering can be observed. The shapes 
of the differential cross sections arc quite different at res- 
onance energies compared to collision energies away from 
the resonances. 

A simple way to analyze the resonances and to quali- 
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in collisions with He, as defined in Eq. ([9]) of the text, for 
total angular momenta J — 1/2 - 11/2. 



tatively explain backward scattering appearing in differ- 
ential cross sections is to calculate the collision lifetime, 
which is the difference between the time that the colli- 
sion partners spend in each other's neig hborhood with 
and without the interaction [jn], [49|, [501 • For a direct 
comparision with the partial cross sections shown in Fig. 
II H we compute the collision lifetime from initial state 7 
to final state 7' for individual total angular momenta J, 
defined as 



Ad y (E)=Re 



-in £ 

l,P,l' .p' 



'VP 



■y',1,1' 



^7,7',!,/' 

dE 



(9) 

where p and I' , p' denotes the orbital angular momen- 
tum and parity of initial and final levels, respectively, 
and 1 1 v denotes S'-matrix elements for total angu- 
lar momentum J from close-coupling calculations. The 
lifetimes vs. J for the j = 3/2, F x , f -> j = 1/2, F 2 ,f 
transition are plotted in Fig. Q2J Clearly, the resonance 
peaks in Fig. [TTJ are well reproduced, with collision life- 
times of a few picoseconds. The most intense resonance 
peak lies at 186.4 cm -1 , which was largely due to the 
,7 = 9/2 partial cross section. From Fig. [T^] we see that 
the corresponding lifetime is about 6 ps. We can com- 
pare this collision lifetime with the rotational period of 
the OH-He van der Waals complex. We estimate the ro- 
tational constant of the complex to be 6.4 x 10 -24 J from 
the expectation value of 1/R 2 computed with the wave 
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FIG. 13: State-to-state inelastic scattering cross sections of 
OH (X 2 H 3 /2, j = 3/2, /) radicals with Ne atoms as a function 
of the collision energy. Cross sections for inelastic collisions 
populating the (a) j = 3/2, Fx, e (black) and the j = 5/2, Fi,e 
(purple) states; (b) the j = 5/2, Fi, f (orange), j = 1/2, F 2 , e 
(brown), j = 3/2, F%, f (pink), and the j — 7/2, Fi, e (green) 
states; (c) the j = 1/2, F 2 , f (red), j = 3/2, F 2 ,e (blue), and 
the j = 7/2, Fx, f (cyan) states. 



function obtained from the DVR method on the lowest 
lying J = 9/2 adiabatic bender curve. This corresponds 
to a rotational period 14.9 ps, assuming / = 3. 

We thus conclude that the collision lifetime has the 
same order of magnitude as the rotational period of the 
OH-Hc complex. It is therefore not surprising to observe 
significant backward scattering at some resonance ener- 
gies. At off-resonance energies the collision lifetime will 
be < 1 ps, which is much smaller than the OH-Hc rota- 
tional period. Hence, backward scattering is expected to 
be barely observable. 



V. OH-NEON COLLISIONS 

To describe the i nteraction betwe en OH and Ne, 
we used the PES of ISumivoshi et all [H. This PES 
was calculated using an explicitly correlated, spin- 
unrestricted coupled-cluster approach [UCCSD(T)-F12b] 
with a quintup le-zeta basis set (AV5Z). Although 
ISumivoshi et al.l calculated a three-dimensional poten- 
tial, we used in Ref. [Hf their interaction potential eval- 
uated at the equilibrium distance r e = 1.832 ao for the 
scattering calculations, so that no effect of the OH vi- 



brational motion was included. In that reference it was 
shown that this procedure already gives excellent agree- 
ment between theory and high-precision scattering exper- 
iments for OH-Ne collisions. Since we found in Section 
IIIII that the vibrational motion of OH can be of quantita- 
tive influence, we also calculated an ad iabatic potential 
from the three-dimensional potential of ISumivoshi etaTl 
in the same way as we did for OH-He. The resulting 
adiabatic potential was found to improve slightly the ex- 
cellent agreement with the experimental results for the 
scattering of OH and Ne. In the present study, we 
use the adiabatic potential throughout and compute the 
cross sections on a much finer grid than in the study by 
IScharfenberg et all (34| in order to study scattering res- 
onances. In Fig. [13] we show the energy dependence of 
state-to-state integral cross sections for collisions of the 
OH radical with Ne atoms, where the OH radicals are 
initially in the j = 3/2, Fx, f level. Overall, the behav- 
ior of the inelastic cross sections as a function of energy 
is rather similar to what was observed for the OH-Hc 
system in Section IIVI For example, we again observe 
a propensity for transitions preserving the total parity. 
However, in the OH-Ne system none of the channels ap- 
pears to have particularly strong Feshbach resonances, 
as was the case for the j = 3/2, jFi, / — > j = 1/2, F2, f 
transition of the OH-He system. The most pronounced 
resonant features observed for OH-Ne collisions are shape 
resonances in the j = 3/2, Fx, f —> j = 5/2, Fi, / transi- 
tion. In Fig. [T3]we show these shape resonances in more 
detail. 

Looking closely at Fig. [H] we see several resonance 
peaks that correspond to an increase in the cross section 
by about a factor of two compared to the nonresonant en- 
ergies. A relatively strong resonance occurs at a collision 
energy of 99.23 cm" 1 ; this resonance increases the cross 
section by a factor of four compared to the background. 
The latter resonance is indicated by the Roman numeral 
V. The main contributions to this resonance originate 
from partial cross sections with total angular momenta 
of J = 37/2 and J = 39/2. In Fig. QU we also show 
the differential cross sections for several energies that are 
marked by Roman numerals in panel (a). For the res- 
onances at collision energies of 86.83, 94.90 and 99.23 
cm -1 , the cross sections are shown in the panels (b), (d) 
and (f). In these plots, large amplitudes for backscat- 
tcring are observed. To compare, the differential cross 
sections were also calculated away from the resonances 
at the energies 93.00, 98.00 and 101.00 cm -1 , and the re- 
sults are shown in the panels (c), (e) and (g). In the case 
of nonresonant scattering, the observed backscattering is 
significantly reduced. The differential cross sections at 
these resonances look similar to those at the shape reso- 
nances for the j = 3/2, Fx, f — > j = 5/2, Fx, e transition 
of the OH-He system (see Fig. [7]), where also an in- 
crease in backscattering was found. With a measurement 
of the differential cross sections, the strong increase and 
decrease in the backscattering might help in experimen- 
tally identifying the shape resonances at 94.90 and 99.23 
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FIG. 14: (a) Integral cross section for collisions of OH radicals with Ne atoms as a function of the collision energy. The initial 
state of the OH radical is the j — 3/2, J 7 !, / state and the final state is the j = 5/2, Fi,/ state. The blue curve shows the 
theoretical results without any convolution, the red curve shows the integral cross section convoluted with a Gaussian energy 
distribution having a FWHM of 1 cm -1 , while for the green curve a FWHM of 5 cm" was used. At the collision energies 
designated with the Roman numerals I- VI, differential cross sections were calculated, shown in panels (b) to (g). Again the 
blue curves are not convoluted, while for the red curves a FWHM of 1 cm -1 was used, and for the green curves a FWHM of 5 
cm -1 At the scattering resonances (I, III and V), strong backscattering is observed. 



cm . However, we note that one must be careful with 
identifying the backscattering signal with resonances in 
the cross section. Namely, in Fig. H4Ta) also less pro- 
nounced resonances are seen, and not all of them have 
such a strong backscattering signal as the strongest reso- 
nances I, III and V. Moreover, closer to the threshold of 
the j = 3/2, Fl, / — > j = 5/2,Fi,e transition significant 
backscattering is observed away from the resonances. 

VI. DISCUSSION 

The experimental observation of resonance structures 
as discussed in this paper would comprise a very detailed 
test for the calculated PES's and scattering calculations 
on these PES's. The Stark deceleration technique pro- 
vides a source of state-selected molecules with a tunable 
velocity and narrow velocity distribution. This technique 
enables state-to-state scattering experiments in which the 
collision energy can be precisely tuned over a wide range 
with a high collision energy resolution. Yet, the obser- 



vation of scattering resonances requires an energy res- 
olution that has not yet been achieved in this type of 
experiments. In this section, we analyze the collision en- 
ergy resolution required to observe resonance features in 
either the statc-to-statc integral or the differential cross 
sections for OH-He and OH-Ne collisions. We discuss 
the requirements on the beam velocity and angular dis- 
tributions, and discuss the feasibility of obtaining these 
distributions. 

Referring back to Figs. [Tl [TQl and[T4l the most promi- 
nent resonance structures are found for OH-He in the 
j = 3/2, j = 3/2, F^e, j = 3/2, Fi,/ -»• 
j = 5/2, F u e and the j = 3/2, Fi,/ -> j = 5/2, F 2 J 
transitions. For OH-Ne collisions the j = 3/2, Ft, / — > 
j = 5/2, Fl,/ channel is of most relevance. To simulate 
what would be observed in a molecular beam scattering 
experiment, the integral cross sections in Figs. [71 QUI and 
1141 are convoluted with Gaussian collision energy distri- 
butions of 1 and 5 cm -1 full width at half maximum 
(FWHM). In Fig. [TSJ the resonance structure at colli- 
sion energies around 85 cm -1 for the j = 3/2, Fl, / — > 



12 



4 









1 1 


5 cm" 1 


' ? cm- 1 






1 cm- 1 






0.5 cm- 1 




aHX/V 


theoretical 
cross section 


i i i i 



70 75 80 85 90 95 

collision energy / cm- 1 



FIG. 15: Inelastic scattering cross section for the j = 
3/2, F%, f — ¥ j = 3/2, F\, e channel in collisions of OH radicals 
with He atoms, showing Feshbach resonances that correspond 
to the opening of the j — 5/2, Fi rotational levels of OH. The 
theoretical curve is convoluted with Gaussian energy distribu- 
tions with FWHM of 0.5 cm" 1 , 1 cm" 1 , 2 cm" 1 , and 5 cm" 1 . 
The upper four curves have been given a vertical offset for 
reasons of clarity. 



j = 3/2, Fi, e channel in OH- He collisions is shown. This 
scattering channel displays a number of Feshbach reso- 
nances, corresponding to the opening of the j = 5/2, Fx 
channels, that are grouped within a relatively narrow 
range of collision energies. The theoretical curve is convo- 
luted using 0.5, 1, 2, and 5 cm" 1 (FWHM) energy distri- 
butions. We will use this scattering channel as a bench- 
mark to establish the energy resolution required in the ex- 
periments to observe signatures of scattering resonances. 
Scattering resonances are partially resolved for energy 
resolutions of < 1 cm" 1 . When a resolution between 1 
and 2 cm" 1 is achieved, some of the resonance structure 
is resolved, while for resolutions of 5 cm" 1 , most of the 
resonance structure has disappeared. In these cases, at 
best only a broad peak is observed in the integral cross 
section. 

For given beam velocity and angular spread, the col- 
lision energy distribution is a function of the collision 
energy; the highest resolutions are obtained for the low- 
est collision energies. We can estimate the beam pa- 
rameters required to reach collision energy distributions 
of < 2 cm" 1 at a collision energy of 85 cm" 1 for col- 
lisions of Stark-decelerated OH radicals with He atoms, 
i.e., to (partially) resolve the scattering resonances shown 
in Fig. \T5\ We assume a beam intersection angle of 
45° and choose the velocities of the He and OH beams 
such that the relative velocity vector is perpendicular to 
the He atom velocity vector. In this geometry, the colli- 
sion energy distribution is almost independent of the He 
atom beam velocity spread (l6| . A collision energy of 
85 cm -1 is reached for He and OH velocities of 790 m/s 
and 1120 m/s, respectively. In the chosen geometry, the 
most critical parameter that determines the collision en- 
ergy resolution is the distribution in beam intersection 
angles Atp. If we assume extremely well collimatcd beams 



such that Acj) = 10 mrad (corresponding to 0.6°), a ve- 
locity spread of the OH radicals of 5 m/s results in a 
collision energy resolution of 1.9 cm" 1 . 

Experimentally, the most challenging requirement is 
the angular spread of both beams. Multiple collimation 
slits for both the OH and He beams are required to reach 
the required angular spreads. The required He atom ve- 
locity can be obtained using a cryogenic source that is 
maintained at a temperature of 60 K, and the required 
OH velocity can be produced using the Stark decelerator. 
The required velocity spread for the OH radicals can be 
obtained using the Stark decelerator, either by choosing 
the appropriate phase angle in the decelerator [Hil l , o r 
by additional phase-space manipulation techniques [52j. 
For the OH-Ne system, even more stringent requirements 
apply to the beam distributions due to the higher re- 
duced mass for this system. Beam speeds of 664 m/s 
and 470 m/s for the OH radical beam and Ne atom, re- 
spectively, will result in a collision energy of 85 cm" 1 
using a beam crossing angle of 45°. For the velocity and 
angular distributions used above for OH-He, a collision 
energy distribution of 2.1 cm" 1 is obtained. 

The signatures of scattering resonances can also be in- 
ferred from differential cross sections. State-to-state dif- 
ferential cross sections can be measured using the veloc- 
ity map imaging technique, that provides the full angu- 
lar and velocity distribution of the scattered molecules 
[53l ]. Alternatively, information on the differential cross 
section may be obtained via Doppler profile measure- 
ments of the scattered molecules. The backward scat- 
tered components that appear in the differential cross 
sections when a resonance is accessed offers interesting 
prospects to reveal the existence of resonances. As the 
collision energy is tuned over a group of scattering res- 
onances, the presence of these resonances can in princi- 
ple be inferred from the measured product flux at back- 
ward scattering angles. In Figs. [TUl and HU the 
differential cross sections are shown at collision ener- 
gies near and at the resonances, convoluted with Gaus- 
sian collision energy distributions of 1 and 5 cm" 1 full 
width at half maximum (FWHM). In particular for the 
j = 3/2, Fi,/ j = 5/2, Fi,/ channel in OH-Ne, a 
significant scattering intensity at backward scattering an- 
gles remains, even for a collision energy resolution as high 
as 5 cm" 1 . In the integral cross section, no signature of 
the scattering resonances is observable at these energy 
resolutions. In these cases, it may be favorable to exper- 
imentally explore the existence of resonances via mea- 
surements of differential cross sections instead of integral 
cross sections. 

The examples treated above show that collision energy 
resolutions in the 1-2 cm" 1 range, although very chal- 
lenging experimentally, should allow for the observation 
of both shape and Feshbach resonances in the integral 
cross sections for inelastic collisions between OH radicals 
and He or Ne atoms. The OH radical is an excellent 
candidate in these experiments, as a sensitive detection 
scheme, appropriate for ion imaging, has recently been 
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developed for this species [54| . In addition, the relatively 
large rotational spacing of the molecular levels results 
in a molecular beam pulse with less initial population in 
excited rotational states, and therefore a packet of Stark- 
decelerated OH radicals with a high state purity. This 
facilitates the sensitive and background free detection of 
scattering products and enables the implementation of 
beam collimators that improve the angular and velocity 
spreads of the beams at the cost of particle densities. The 
disadvantage of the large rotational spacing, however, is 
the relatively high energies of the energetic thresholds 
for inelastic scattering, and corresponding relatively high 
collision energies at which scattering resonances appear. 
In this respect, the inelastic scattering of OD radicals 
or ND3 molecules with He and Ne atoms will be inter- 
esting candidates for studying scattering resonances as 
well, reducing the energy for the lowest lying threshold 
to 42 cm -1 and 14 cm -1 , respectively. To reach a colli- 
sion energy resolution < 2 cm -1 at these energies will re- 
lax the requirements on the velocity and angular spreads 
of both beams. We will investigate resonance effects in 
the rotationally inelastic scattering of OD and ammonia 
molecules with He atoms in forthcoming publications. 

VII. CONCLUSIONS 

We have presented detailed calculations on scattering 
resonances in the rotationally inelastic scattering of OH 
radicals with He and Ne atoms. For OH-He, we have de- 
veloped new 3D ab initio potential energy surfaces, and 
the inelastic scattering cross sections that are derived 
from these surfaces compare favorably with recent ex- 
periments. We have identified numerous scattering reso- 
nances - of both the shape and Feshbach types - in the in- 
tegral cross sections. We have analyzed these resonances 
using the adiabatic bender model and computed collision 
lifetimes. We observe dramatic changes in the differen- 



tial cross sections at the resonances, showing in selected 
cases a forward-backward peaking of the scattered flux. 
The analysis of scattering resonances presented here will 
be indispensable in the experimental search for such res- 
onances in, for instance, crossed beam scattering experi- 
ments using Stark-decelerated molecular beams. To ex- 
perimentally observe signatures of resonances in the in- 
tegral cross sections and to partially resolve individual 
resonances, a collision energy resolution of < 2 cm -1 is 
required. Obtaining energy resolutions < 2 enr 1 mainly 
requires highly collimated molecular beams, which ap- 
pears challenging. Alternatively, signatures of scattering 
resonances may be found in the differential cross sections. 
The selective detection of scattered molecules at back- 
ward scattering angles may facilitate the identification of 
resonances if the collision energy resolution is not suffi- 
cient to resolve them in the integral cross sections. 
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